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ABSTRACT 


Pulse-Doppler  processing  is  an  often  used  technique  for  extracting  range 
and  Doppler  information  of  targets.  Such  a  system  can  be  described  as  con¬ 
sisting  of  three  basic  components,  a  pulse  compressor,  sampler,  and  Doppler 
processor.  When  quadrature  detection  is  used,  the  gains  and  timing  of  the 
sampled  lnphase  and  quadrature  components  of  the  pulse  compressor  output  must 
be  matched  in  order  to  minimize  distortion  and  noise  at  the  Doppler  processor 


output.  In  addition,  the  jitter  in  the  sampling  times  also  causes  distortion 
and  noise.  This  report  presents  an  analysis  of  these  effects  and  describes 


the  signal  degradations  due  to  mismatch  and  jitter  in  the  sampler. 
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1.  INTRODUCTION 


Pulse-Doppler  processing  is  used  to  obtain  range  and  velocity  informa¬ 
tion  of  targets  for  radar  applications  such  as  ballistic  missile  defense.  A 
radar  system  based  on  this  type  of  processing  transmits  a  waveform  which  con¬ 
sists  of  a  sequence  of  pulses  uniformly  spaced  in  time.  Each  of  these  pulses 
is  linear  frequency  modulated  (LFM)  and  is  identical  to  the  other  pulses  (ex¬ 
cept  possibly  for  a  weighting  factor).  The  components  of  this  waveform  are 
shown  in  Figure  1. 

The  received  waveform  (for  a  single  target)  is  very  similar  to  that 
shown  in  Figure  1,  except  that  the  signal  is  delayed  in  time  corresponding  to 
the  target  range,  and  the  signal  is  Doppler  Bhifted  in  frequency  corres¬ 
ponding  to  target  velocity.  In  addition,  noise  is  added  to  the  signal.  The 
received  signal  is  first  filtered  by  a  pulse  compressor  (matched  filter  for  a 
LFM  pulse),  producing  the  signal  shown  in  Figure  2.  This  pulse  compressed 
signal  is  converted  from  a  continuous  signal  to  a  sampled  data  set  by  a  sam¬ 
pler,  and  the  sampled  data  is  then  processed  by  a  Doppler  processor.  For  a 
given  range  cell,  the  Doppler  processor  computes  the  discrete  Fourier  trans¬ 
form  (DFT)  of  a  set  of  samples  rtiich  are  spaced  in  time  by  an  amount  equal  to 
the  pulse  spacing  shown  in  Figures  1  and  2.  The  number  of  points  used  in  the 
DFT  is  equal  to  the  number  of  pulses  in  the  waveform.  The  Doppler  shift  of 
the  target  is  proportional  to  the  rate  of  change  of  phase  from  sample  to 
sample,  and  appears  in  the  DFT  output  as  a  peak  at  the  Doppler  frequency. 
Figure  3  shows  a  simplified  diagram  of  a  pulse-Doppler  processor. 
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Pigur*  2.  Pulse  Compressed  Waveform 


Figure  3.  Pulse-Doppler  Processor 
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One  of  the  limitations  in  the  system  just  described  is  that  the  sampler 
is  not  perfect,  but  instead  exhibits  a  certain  amount  of  mismatch  between  the 
gains  and  sampling  times  of  the  inphase  and  quadrature  (I  and  Q,  real  and 
imaginary)  components  of  the  sampled  pulse  compressed  signal.  The  sampling 
times  are  also  subject  to  random  jitter.  These  nonideal  characteristics 
cause  noise  and  distortion  to  appear  in  the  DFT  output,  so  that  target 
detection  becomes  more  difficult. 

The  purpose  of  this  paper  is  to  describe  analytically  the  degradation  of 
pulse-Doppler  processing  due  to  sampler  Imperfections.  This  analysis  re¬ 
quires  an  understanding  of  the  signal  being  sampled,  the  pulse  compressed 
waveform  shown  in  Figure  2.  This  waveform  can  be  described  in  terms  of  the 
time-frequency  autocorrelation  function. 

The  time-frequency  autocorrelation  x(t»f)  of  a  signal  p(t)  is  obtained 
by  convolving  a  frequency  shifted  and  time  shifted  version  of  p(t)  with  a 
matched  filter  for  p(t).  From  Nathanson  [1,  p.  284], 

00 

x(t,f)  »  /  p(u)p*(u+t)e  ^2lTfudu  . 

— « 

For  example,  the  formula  of  a  LFM  pulse  [1,  pp.  287-288]  is 

p(t)  -  explj2rr(f0t  +  kt2/2)],  jtj  <  x/2  (1) 

for  pulse  length  t  and  bandwidth  kx ,  and  the  corresponding  function  x(t,f)  is 

jitft  sln[n(kt+f ) (x—  [ t j > ]  ,  jtj  <  x  .  (2) 

X<t’f>  “  6  1T(kt+f)T  '  '  M  ” 

Note  that  in  Equation  (2),  x(t.f)  "  x(-t*“f)*  The  nature  of  x(i»i)  is  that 
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there  Is  a  peak  at  t  ■  -f/k  with  nominal  width  1/kx.  Away  from  the  peak, 
x(t,f)  diminishes,  and  Is  Identically  zero  for  |y|  >  t.  The  purpose  of  pulse 
compression  Is  to  obtain  a  width  1/lcr  which  Is  much  less  than  x.  The 
narrowing  of  the  width  yields  Improved  range  resolution. 

The  transmitted  radar  waveform  is  a  sequence  s(t)  of  LFM  pulses 

N-l 

s(t)  -  l  p(t-mT)  (3) 

m*0 

where  T>t  Is  the  pulse  spacing  indicated  in  Figure  1.  The  received  radar 
signal  r(t)  Is  a  Doppler-shifted  and  time-delayed  version  of  s(t). 

r(t)  -  s(t-p)e^2wft 

Letting  p*(-t)  be  the  pulse  compressor,  the  pulse  compressed  waveform  is 
obtained  by  filtering  r(t)  with  p*(-t).  The  result  is 

N-l 

y(t)  -  l  x(t-mT-p ,f )e^  mT  +  n(t)  (4) 

m«0 

where  n(t)  is  a  random  noise  signal,  p  Is  the  delay  due  to  range,  and  f  is 
the  Doppler  shift.  For  simplicity,  we  are  assuming  a  target  signal  amplitude 
of  unity.  In  a  pulse-Doppler  radar,  the  width  of  the  peak  In  x(t,f)  Is  much 
less  than  t  and  thus  much  less  than  T.  Thus,  we  can  make  an  approximation 
for  y(t)  In  the  vicinity  of  the  peaks, 

y(t)  -  x(t-mT-p ,f +  n(t)  (5) 

for 

(mT+p-f/k)  -  l/2kx  <  t  <  (mTfp-f/k)  +  l/2kx 
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This  approximation  is  actually  an  equality  when  T  >  T+l/2kx  since  x(t*f)  is 
zero  for  |t|  >  t.  For  an  N  pulse  waveform,  Doppler  processing  is  performed 
by  computing  the  DFT  of  N  samples  of  y(t)  uniformly  spaced  with  spacing  T. 

If  we  let  the  m-th  sample  occur  at  t  ■  t'-f/k+mT+p  where 
j  t *  |  <  l/2fcr ,  then  the  DFT  input  samples  are 

xm  “  x(t’-f/k,f)ej2llf,DT  +  n(tm)  ,  m=0 . N-l  (6) 

and  each  of  the  samples  occurs  within  the  peak  of  x»  This  •  mition  of  t' 
corresponds  to  examining  the  range  cell  in  which  the  target  where  t'  is 

the  distance  from  the  center  of  the  range  cell.  The  analys  ~  the  fol¬ 
lowing  sections  is  made  easier  by  writing  x  in  terms  of  the  .^al-valued 

in 

function  G(t,f ) , 


G(t,f)  -  x(t,f)e-jlTft 


sin[ir  1 

Ikt+f )(t- 

[>1 

IT  < 

!kt+f )t 

1 

From  Equations  (6)  and  (7),  we  have 

x  =  G(t’-f/k,f  )ej2l'fmT+j0  +  n(t  ) 
m  m 


(7) 


(8) 


where  0  *  nf(t’-f/k)  and  does  not  depend  on  m.  Note  that  the  samples  in 
Equation  (8)  are  of  a  noisy  sinusoid,  where  the  sinusoid  component  is 
multiplied  by  a  term  which  is  independent  of  the  sample  index  m.  Thus,  the 
DFT  of  xQ,...,xN_1  is  the  DFT  of  a  noisy  sinusoid.  The  location  of  the  peak 
in  the  DFT  output  depends  on  f  and  T,  and  indicates  the  amount  of  Doppler 
shift. 
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In  the  following  sections.  Equation  (8)  will  be  used  to  determine  the 


effects  of  sampling  gain  mismatch,  timing  mismatch,  and  timing  jitter. 
Section  2  investigates  gain  mismatch,  Section  3  examines  sample  time 
mismatch,  and  Section  4  discusses  the  effects  of  sample  timing  jitter. 


2.  GAIN  MISMATCH 


This  section  examines  the  effect  of  mismatched  I  and  Q  gains  on  the  sam¬ 
pling  process.  It  will  be  shown  that  the  result  is  the  appearance  of  a  false 
target ,  plus  reduced  SNR. 

If  we  let  a  and  b  be  the  gains  of  the  I  and  Q  components  of  the  sampler, 
then  the  samples  of  Equation  (8)  become 

»  [a  cos(2irfmT+0)  +  jb  sin(2irfmT+Q) )G( t*  -f /k,f  ) 

+  an^m)  +  jbnQ(m)  (g) 

where  n^Xm)  and  n^m)  are  the  I  and  0  components  of  the  noise.  The  form  of 

can  be  arranged  to  show  the  resulting  degradations  without  having  to 

compute  the  DFT  of  these  samples.  Equation  (10)  shows  that  x  is  equivalent 

m 

to  samples  of  a  signal  which  consists  of  two  targets  plus  noise,  where  the 
targets  are  at  the  same  range,  but  with  opposite  Doppler  shifts. 

x  -  1&.  eJ2"fmT+J0  +  l^b  e-j2w£mT-j0jG(t,_f/k  f) 
tn  l  t 

+  anj.  (m)  +  jbn^Cm)  (10) 

The  first  term  in  Equation  (10)  corresponds  to  t  te  true  target,  and  the 
second  term  is  a  false  target.  The  ratio  of  the  false  target  signal  power  to 
the  true  target  signal  power  (FTR)  is  easily  computed  as  the  ratio  of  the 
square  of  the  relative  amplitudes. 
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(11) 


The  relative  strength  of  the  false  target  Is  shown  in  Figure  4  as  a  function 
of  the  ratio  of  I  and  Q  gains. 

The  loss  In  slgnal-to-nolse  ratio  (SNR)  due  to  gain  mismatch  can  also  be 
easily  determined.  If  the  noise  samples  n^(m)  and  n^m)  are  Independent  and 
have  equal  variances  (power),  then  the  noise  power  Is  multiplied  by  the 
factor 


2  2 
a  +b 


(12) 


where  g  is  the  Ideal  gain.  Similarly,  the  signal  power  Is  multiplied  by  the 
factor 


(a+b)' 

**2 


(13) 


The  ratio  of  these  two  functions  gives  the  ratio  between  the  mismatched  SNR 
and  the  matched  SNR,  which  Is  a  net  loss.  This  loss  In  db  Is 


LOSS  - 


10 


log 


2(a2+b2) 

(a+b)2 


(14) 


Figure  5  shows  the  loss  as  a  function  of  gain  mismatch.  It  Is  seen  that 
the  loss  is  not  large  for  reasonable  amounts  of  mismatch.  The  next  section 
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R  loss 


discusses  the  effects  of  timing  mismatch.  It  will  be  shown  that  timing 
mismatch  can  be  related  to  a  target  dependent  form  of  gain  mismatch,  again 
resulting  in  a  false  target  and  SNR  loss. 
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3.  SAMPLING  TIME  MISMATCH 

In  this  section,  we  investigate  the  result  of  timing  mismatch  between 
the  I  and  Q  samples.  As  In  the  previous  section,  the  analysis  will  be  re¬ 
stricted  to  a  single  target  In  the  presence  of  noise.  The  results  carry  over 
to  multiple  targets. 

To  investigate  the  effects  of  timing  mismatch,  we  modify  Equation  (8)  so 

that  the  I  component  is  sampled  at  time  tm+<S  and  the  Q  component  is  sampled 

at  time  t  -6  where 
m 

t  -  t'-f/k+mT+p  (15) 

m 

is  the  nominal  sampling  time  for  the  rn-th  sample.  In  this  case,  the  samples 
are 

x  “  cos[2nf  (mT+-5)-K5]G(t’-f/k+6  ,f  ) 
m 

+  jsin[2ir  f  (raT-6  )+6  JG(  t’  -f /k-6  ,f  ) 

+  n;[(m)  +  jnq(m)  (16) 

where  n^  and  n^  are  the  time-shifted  noise  terms.  If  the  two  noise  terms  are 
independent  and  represent  a  stationary  random  process,  then  the  noise  power 
and  statistics  are  independent  of  6 .  In  Appendix  A,  it  is  shown  that  the 
effect  of  6  in  the  sin  and  cos  functions  of  Equation  (16)  is  much  less 
significant  than  that  of  6  in  the  function  G,  and  that  the  6  in  the  sin  and 
cos  can  be  dropped. 
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As  In  Section  2,  the  expression  for  the  samples  can  be  arranged  to 


appear  as  two  targets  plus  noise: 


.  j2itfmT+j0 

x  -  4  [G(t’-f/k+6 ,f)  +  G(t'-f/k-6 ,f )]e 
m  z 

,  -j2nfraT-j0 

+  j  [GCt'-f/k-hS  ,f)  -  G(t'-f/k-6,f  )]e 


+  n^m)  +  jnQ(m) 


(17) 


The  first  term  corresponds  to  the  true  target  and  the  second  term  corresponds 
to  a  false  target  at  the  same  range,  but  with  the  opposite  Doppler  shift. 

The  ratio  of  false  target  power  to  true  target  power  is 


m  „  t G( t * -f /k+6  ,f  )  -  G(t'-f/k-6,f)]2 
lG(t'-f /k+6  ,f  )  +  G(t’-f/k-6,f  )]2 


(18) 


Since  the  noise  power  is  unaffected  by  timing  mismatch,  the  loss  in  SNR  is 
equal  to  the  loss  in  signal  power.  In  db  relative  to  the  no  mismatch  case, 
this  loss  is 


LOSS  -  10  log 


4G(t’-f/k,f)‘ 


(G(t’-f/k+6  ,f)  +  G(t’-f/k-6 ,f )] 


(19) 


Notice  that  both  the  FTR  and  the  SNR  loss  depend  strongly  on  the  target  range 
and  Doppler  shift.  Figure  6  shows  the  FTR  as  a  function  of  the  mismatch  26 
(expressed  as  a  fraction  of  the  Nyquist  sampling  interval  1/kt)  for  the  case 
where  the  bandwidth,  kr ,  is  much  greater  than  the  Doppler  shift,  f,  and  where 
the  pulse  length,  t,  is  much  longer  than  the  sampling  interval,  1/kx.  The 
worst  case  is  obtained  when  the  peak  of  the  response  is  located  between  two 
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range  cells  (l.e.,  the  peak  Is  located  1/2  sampling  period  from  where  the 
sample  Is  actually  taken).  With  Nyqulst  sampling,  the  sampling  period  is 
1/kx  and  the  worst  case  occurs  when  t'  ■  ±  l/2kx.  The  derivation  is  not 
given  here,  but  it  can  be  shown  that  the  worst  case  is  approximately 
20  log(26kx).  The  average  case  is  obtained  by  averaging  the  FTR  for  target 
peaks  located  from  1/2  sampling  period  before  the  nominal  sampling  point  to 
1/2  period  after.  It  is  interesting  to  note  that  the  average  is  a  constant  6 
db  better  than  the  worst  case.  Figure  7  shows  the  SNR  loss  for  the  example 
of  Figure  6. 

The  above  analysis  can  easily  be  adapted  to  include  both  gain  and  timing 
mismatch.  In  this  case,  we  have 


LOSS  -  10  log 


2(a2+b2)G(t,-f/k.f )2 _ 

[aG(t'-f /k+6  ,f  )  +  bG(t'-f/k-6,f  )]‘ 


Thus  far,  we  have  seen  that  both  gain  mismatch  and  timing  mismatch  lead 
to  false  targets  and  reduced  SNR.  The  next  section  examines  the  effects  due 
to  timing  jitter.  It  will  be  seen  that  jitter  causes  additional  target 
dependent  distortion. 
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Figure  6.  False  Target  Level  for  Timing  Mismatch 
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Figure  7.  SNR  Loss  for  Timing  Mismatch 


4.  TIMING  JITTER 


In  this  section,  we  examine  the  effects  due  to  sampling  time  jitter.  Two 
different  models  are  considered  for  jitter.  The  case  of  frequency  jitter  In 
the  system  sampling  clock  Is  discussed  first,  and  Is  followed  by  a  discussion 
of  Independent  random  jitter.  The  final  part  of  this  section  combines  both 
jitter  models. 

4.1  Clock  Frequency  Jitter 

The  samples  In  the  case  of  sampling  frequency  jitter  are  modelled  as 

sampling  the  I  channel  at  time  u  «t  46+r  with  gain  a,  and  the  Q  channel 

m  m  m 

at  time  v  «t  -6+r  with  gain  b.  The  terms  a,b,6  represent  the  gain  and 
m  in  m 

timing  mismatch  discussed  In  earlier  sections.  The  r  are  a  random  process 

m 

representing  the  jitter. 

A  good  model  of  this  jitter  process  is  to  represent  the  difference 

2 

rm"rm-l  as  an  independent  Gaussian  process  with  zero  mean  and  variance  o  . 

The  justification  for  this  model  Is  that  the  system  clock,  which  Is  used  to 
drive  the  sampling  process,  can  be  represented  as  having  a  random  clock 
period  which  is  independent  from  one  clock  cycle  to  the  next,  and  has  a 
stationary  mean.  The  system  clock  operates  at  a  frequency  greater  than  or 
equal  to  the  sampling  rate,  which  in  turn  is  greater  than  or  equal  to  the 
bandwidth,  kx  .  The  samples  used  for  the  Doppler  processing  of  a  given  range 
cell,  however,  are  spaced  by  the  pulse  spacing,  T,  as  described  In  Section  1. 
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Since  the  pulse  spacing  is  much  greater  than  the  sampling  period,  the  pulse 
spacing  is  the  sum  of  a  large  number  of  system  clock  periods.  Since  these 
clock  periods  are  modelled  as  independent,  identically  distributed  random 
variables,  we  have  by  the  central  limit  theorem  [21  that  and 

vm-vm_i,  the  time  between  two  samples  used  for  Doppler  processing,  are 
Gaussian  distributed  random  variables,  with  mean  T.  (The  variance  of  these 
random  variables  is  defined  as  o^.  Note  that  is  the  product  of  the 
variance  of  the  system  clock  period  and  the  number  of  clock  periods  which 
occur  during  time  T.  Thus,  is  proportional  to  T. )  By  the  definition 
given  above  for  and  vm,  we  have  that  rm-rm_j  is  Gaussian  with  zero 
mean  and  variance  o^. 

As  an  example,  consider  the  case  of  kr“50  MHz  bandwidth,  and  T-4  psec 
pulse  spacing.  The  number  of  system  clock  periods  between  two  samples  used 
for  Doppler  processing  is  at  least  krT-200. 

Without  loss  of  generality,  we  define  p,  the  time  delay  due  to  target 
range,  such  that  the  jitter  is  zero  at  time  (to+tfj_i)/2.  (As  in  the 
previous  sections,  only  the  range  cell  which  includes  the  target  is 
analyzed.)  Since  Unf-u,,,-!  and  v^v^i  are  Gaussian  with  mean  T  and 
variance  o^,  the  defintion  of  p  implies  that  ^^-(uo+Ufj-i )/2  and 
Vm-(V0+VN-1 )/2  are  Gaussian  with  mean  mT-T(N-l)/2  and  variance 
|m-(N-l)/2 |o^.  Thus,  rm  is  Gaussian  with  zero  mean  and  variance 
|m-(N-l)/2|o2. 

In  order  for  the  Doppler  processing  to  be  effective,  it  is  necessary 
that  the  standard  deviation  of  the  jitter  be  less  than  the  sampling 
interval.  That  is,  o/ ( N— 1 )/2  <  1/kt  for  the  extremes  m*0  and  m-N-l. 


Applying  the  above  definitions  to  Equation  (8)  gives  the  following 
samples, 

x  =  aG(t'-f /k+r  +6  ,f  )cos[2irfmT  +  0+irf(r+6)] 
mm  m 

+  jbG(t'-f/k+r  -6  ,f  )sin[2itfmT  +0  +  irf(r  -6)] 
m  m 

+  an^m)  +  jbnQ(m) 


where  0  =  irf(t'-f/k)  and  where  n^m)  and  n^(m)  are  the  noise  components. 
Similar  to  the  approximation  made  In  the  previous  section,  It  Is  assumed  that 
the  sampling  mismatch  plus  jitter  Is  less  than  one  sampling  period. 

Since  the  sampling  period  Is  no  greater  than  1/kT  ,  we  have  |rm±6 | <l/kr . 

Again,  using  the  result  given  in  Appendix  A,  we  have  that  the  arguments  to 
the  sine  and  cosine  In  Equation  (22)  are  essentially  2TrfraT+0.  The  samples 
can  now  be  written  as 

x  -  aG(t'-f/k+r  +d  ,f  )cos(2irf  mT+0)  +  jbG(t’-f/k+r  -6  ,f  )  sln(2irfmW3) 
m  m  m 


+  anjCnO+jbn^ra) 


(23) 


These  samples  can  be  rewritten  as  the  sum  of  a  jitter-free  signal  and  a 
jitter-dependent  signal. 

xm  -  aG(t'-f/k-W  ,f  )cos(2trfmT+0)  +  jbG(t’ -f /k-6  ,f  )sln(2irfmT+0) 

+  act  cos(2vfmT+0)  +  jb8  sin(2xfmT+0) 
m  m 

+  an  (m)  +  jbn  (m) 

1  0 
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where  the  terms  containing  jitter  are 


a  -  G(t'-f/k+r  -hS  ,f  )-G(t' -f /k+6  ,f ) 
m  m 

B  =  G(t'-f/k+r  -6  ,f)-G(t'-f /k-6  ,f)  .  (25) 

m  m 

From  the  analysis  of  the  preceding  section,  we  know  that  the  jitter-free 
terms  of  Equation  (24)  lead  to  a  false  target  and  reduced  SNR  due  to  the  gain 
and  timing  mismatch.  The  jitter-dependent  terms  in  Equation  (24)  represent  a 
random  distortion.  This  distortion  may  be  computed  in  order  to  determine  the 
signal-to-distortion  ratio  (SDR)  due  to  jitter. 

The  expected  value  of  the  distortion  in  the  samples  Xq,...,x^^  is  de¬ 
fined  to  be 

N-l  „ 

D  -  El  T  laa  cos(2trfmT+G)+jb0  sin( 2irfmT+«)  \  ]  .  (26) 

1  m=0 1  m  ”  1 


If  we  define  the  second  moments  of  a  and  B  as 

m  m 


A 


m 


Et«2  i 

m 


B 

m 


E[B2  ] 

m 


then  the  distortion,  D^,  can  be  written  as 


(27) 


N~1 

T  [a2A  cos2(2xfmT+e)  +  b2B  sin2(2nf mT+O) ) 
®  m 


(28) 
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Unfortunately,  a  closed-form  solution  for  A  and  B  Is  not 

mm 


In  order  to  provide  an  understandable  expression  for  D, ,  A  and 

i  in 

computed  by  approximating  and  Bm  with  a  second-order  Taylor 
the  point  r“0  [3].  Making  .he  definitions 


gl  “  "Hr  G(t,_f/kfr+,s  »£>  r<o 


R2  =  J - 2  G<t,-f/k+r+lS  *£)  |r=0 

dr 

hl  “  h  G<  t'  -f /k+r-6  ,f  )  |  r-0 


*2  "I  T7  G<  t'  ~f  /k+r-6  ,f  )  | 
dr 


gives  the  approximations 


“m  *  Vm  +  *2rra 


B  =  h.r  +  h„r 
ra  Ira  2  m 


available. 

B  may  be 
m 

series  about 


(29) 


(30) 
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Note  that  g^  and  *  and  thus  ,  depend  on  target  Doppler  shift  and  range. 

When  mismatch  Is  also  present,  the  distortion,  ,  can  be  bounded  by 
2  2 

noting  that  cos  (tKl  and  sin  (t)£l.  The  result  Is  a  formula  similar  to 
Equation  ( 34) , 
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(35) 


„  .  I  4,  2  2  2,  2.  ,„3  .  1  2,  2  2^.  2,  2.„2 

D1  <  o  (a  «2+b  h2)(N  -N)  +  -j-  a  (a  g ^+b  n ^ ) N 

In  the  case  of  no  mismatch,  the  upperbound  In  Equation  (35)  Is  two  times  the 
formula  In  Equation  (34). 

At  this  point  all  that  remains  in  computing  the  SDR  is  to  compute  the 
signal  power.  Applying  the  discussion  in  Section  3  to  the  samples  in 
Equation  (24),  the  amplitude  of  the  true  target  is 

j  laG(t’-f/k-H$  ,f )  +  bG(t’-f/k-6,f  ))  .  (36) 

The  signal  power  in  samples  Xq,...,x^_^  is  then 

N-i  ,  2 

S  “It  laG(t'-f/k+6 ,f)  +  bG(t'-f/k-6 ,f )] 
mO 

*  |  (aG(t'-f/k-hS  ,f)  +  bG(t’-f/k-6  ,f))2  .  (37) 

In  the  case  of  no  gain  or  time  delay  mismatch,  the  signal  power  is 

S  -  N  a2G2(t’-f /k,f )  .  (38) 

Defining  SDRj  as  S/Dj  and  using  the  above  equations  for  S  and  Dj ,  it  is 

apparent  that  for  a  fixed  pulse  spacing,  the  SDR^  prior  to  Doppler  processing 

is  inversely  proportional  to  N2,  the  square  of  the  number  of  pulses  in  the 

2 

burst.  Since,  as  mentioned  earlier,  o  is  proportional  to  the  pulse  spacing, 
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2 

T,  the  SDR  is  Inversely  proportional  to  T  also.  These  comments  may  be 

combined  to  show  that  the  SDR^  is  inversely  proportional  to  the  square  of  the 

2  2 

length  of  the  pulse  compressed  waveform,  (N-l)  T  . 

As  an  example  of  the  SDR,  consider  the  case  where  the  bandwidth,  kr  ,  ia 
much  greater  than  the  Doppler  shift,  f,  and  the  pulse  length,  t,  is  much 
greater  than  the  sampling  period,  1/kx .  Figure  8  shows  the  SDR  (averaged 
over  target  range)  defined  by  Equations  (34)  and  (38).  This  curve  is  shown 
as  a  function  of  N,  the  number  of  pulses  in  the  waveform,  and  as  a  function 
of  the  amount  of  jitter,  o,  for  a  N=16  pulse  waveform. 

4.2  Independent  Random  Jitter 

The  samples  in  the  case  of  independent  random  jitter  are  obtained  by 

sampling  the  I  channel  at  time  t^+6+r^  with  gain  a,  and  the  Q  channel  at 

time  t  -6+s  with  gain  b.  As  in  earlier  discussions,  a,b,6  represent  gain 
m  m 

and  time  delay  mismatch.  The  and  are  independent  and  identically 
distributed  random  variables  representing  the  timing  jitter,  and  are  modelled 
as  having  zero  mean  and  variance 

The  justification  of  this  independent  random  jitter  model  is  that  jitter 
can  be  observed  in  electronic  systems,  even  though  the  system  clock  is  stable 
during  the  observation  time.  For  example,  in  digital  systems,  noise  (e.g., 
crosstalk  or  thermal  noise)  can  combine  with  the  nonzero  rise  and  fall  times 
of  logic  signals  to  cause  logic  circuits  to  change  states  at  a  somewhat 
random  time  with  respect  to  the  ideal  time  of  the  state  change. 
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(qp)  a 


The  analysis  of  this  random  jitter  case  parallels  the  analysis  in 
Section  4.1,  producing  a  false  target,  reduced  SNR,  and  target  dependent 


distortion.  The  result  similar  to  Equation  (28)  is  that  the  distortion, 

°2  *  is 

N— 1 

D  =  T  [a2cos2(2irfmT+<E>)E[a2]+b2sin2(2irfmT+e)E[32]]  (39) 

2  **_  m  m 

m=»0 

where 

a  =  G(t’-f/k+r  -h5,f)  -  G( t’ -f /k+6  ,f ) 
m  m 

=  G(t’-f/k+sm-6 ,f )  -  G(t'-f/k-6 ,f)  . 

Using  a  second-order  Taylor  series  as  before,  and  modelling  r  and  8 

mm 

as  independent,  identically  distributed  random  variables  with  zero  mean  and  a 
symmetric  density  function  gives  the  approximations 


E(a2! 

m 

,  2  4  .  2  2 

“  3g2y  +  gjY 

ElB2] 

m 

_,2  4,22 
=  3h2r  +  hxr 

(41) 

where  gj ,  g.,  »  >  h2  are  defined  Equation  (29).  Applying  Equation  (41)  to 
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Equation  (39)  and  assuming  no  gain  or  delay  mismatch  gives 

D2  *  a2N(3g^4+g2Y2)  .  (42) 

When  mismatch  is  present,  a  bound  on  the  distortion  can  be  given  as 

„  ,  2,t/„  2  4,2  2*,  .2„,,.  2  4,,  2  2*  //os 

D2  <  a  N(3g2Y  +g^  )  +  b  N(3h2Y  +h^  )  .  (43) 

The  signal-to-distortion  ratio,  SDR2  -  S/D2,  is  computed  by  using  the 

formula  for  S  in  Equation  (37)  or  (38)  with  Equation  (42)  or  (43).  Since 

both  S  and  D2  are  proportional  to  N,  the  ratio  SDR2  does  not  depend  on  N. 

Note  also  that  SDR2  does  not  depend  on  the  pulse  spacing,  T.  However,  SDR2 

does  depend  on  target  range  and  Doppler  shift. 

As  an  example,  the  SDR2  is  plotted  as  a  function  of  y  In  Figure  9  for 

the  same  example  used  in  Section  4.1  with  Figure  8,  assuming  no  gain  or 

mismatch.  As  before,  SDR2  is  averaged  over  target  range. 

In  the  case  of  independent  random  jitter,  it  can  be  shown  that  the 

distortion  component  in  the  Doppler  processor  output  can  be  described  as 

having  two  parts.  The  first  part  is  related  to  the  means  of  a  and  8  in 

in  m 

Equation  (40),  and  produces  a  nonrandom  distortion  which  has  a  peak  at  the 
Doppler  frequencies  of  the  target  and  false  target.  To  show  this  effect, 


0.0005 


let  a  and  0  be  the  means  of  a  and  0  defined  In  Equation  40.  Then,  from 

m  m 

Equation  (24),  the  jitter-dependent  portion  of  the  samples  which  results  from 

the  means  of  a  and  3  Is 
m  m 

aa  cos(2irfmT+0)  +  jb0  sin(2irfmT+0)  (44) 

which  appears  as  samples  of  sinusoids  at  the  Doppler  shifts  of  the  true 

target  and  the  false  target.  In  the  absence  of  gain  and  delay  mismatch,  we 

have  aa  =  b0  and  Equation  (44)  only  contributes  at  the  Doppler  shift  of  the 

true  target.  Thus,  the  distortion  Is  not  spectrally  flat. 

The  effect  of  this  first  part  Is  shown  In  Figure  10  which  shows  the 

distortion  power  at  the  target  Doppler  shift  relative  to  the  corresponding 

signal  power  as  a  function  of  y*  In  this  example,  we  assume  a  16  pulse 

waveform,  bandwidth  kx  much  greater  than  Doppler  shift  f,  pulse  length  x  much 

greater  than  the  sampling  period  1/kx ,  and  no  gain  or  delay  mismatch. 

It  Is  also  assumed  that  that  r  and  s  are  Gaussian  with  zero  mean  and  varl- 

m  m 

ance  y  *• . 

The  second  part  of  the  distortion  In  the  Doppler  processor  output  appears 
as  zero  mean  complex  random  variables.  Subtracting  Equation  (44)  from  the 
samples  in  Equation  (24)  gives  the  following  jitter-dependent  component  of 
the  samples, 

a(a  -a)cos(2irfmT+G)  +  jb(0  -0  )sln(2nfmTf0)  .  (45) 

m  m 

The  details  are  not  presented  here,  but  it  can  be  shown  that  the  DFT  of 
Equation  (45)  leads  to  zero  mean  complex  Gaussian  random  variables,  with 
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unequal  variances  for  the  real  and  Imaginary  parts.  The  derivation  relies  on 
the  independence  of  the  r^  and  sm,  and  applies  the  Lindeberg  condition  [2] 
of  the  central  limit  theorem  in  evaluating  the  real  and  imaginary  components 
of  the  DFT  outputs.  This  result  is  valid  even  when  the  DFT  inputs  are 
weighted  to  reduce  Doppler  sidelobes. 

4.3  Combined  Jitter  Model 

The  choice  of  which  jitter  model  to  use  depends  strongly  on  radar  system 

parameters  (e.g.,  bandwidth,  pulse  spacing,  number  of  pulses)  and  on 

2  2 

hardware  implementation  parameters,  Y  and  a  .  When  the  product  NT  is  large, 
the  clock  frequency  stability  approach  (Section  4.1)  becomes  more  important 
than  the  independent  random  jitter  approach  (Section  4.2).  However,  when  NT 
is  small,  clock  frequency  stability  is  less  important  than  the  independent 
random  jitter  which  is  always  present. 

In  those  cases  where  neither  jitter  model  is  clearly  dominant,  the  I  and 

Q  timing  jitter,  rm  and  sffl,  may  be  described  as  being  Gaussian  with  zero  mean 

p  ,  i  o 

and  variance  v  '•>  m-(N-l)/2  a  .  This  probability  density  may  be  used  in 
computing  xhe  second  moments  of  offl  and  (5^  in  Equation  (27)  which  are  then 
used  to  compute  the  distortion.  The  results  lead  to  a  distortion  formula  in 
a  straightforward  manner;  the  derivation  is  not  presented  here.  The  result 
is  interesting,  in  that  the  combined  distortion,  D,  is  the  sum  of  and  D2 
plus  a  cross  term.  In  the  case  of  no  gain  or  delay  mismatch, 

D  *  Dj  +  D2  +  |  a2g2Y2o2N2  (46) 
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where  the  Taylor  series  approach  is  used  again,  and  Dj  and  D2  are  described 
by  Equations  (34)  and  (42).  When  mismatch  is  present,  we  have  the  bound 

D  <  B:  +  B2  +  Y2o2N2(a2g2-H)2gi)  (47) 

where  B^  and  B2  are  the  bounds  defined  for  and  D2  in  Equations  (35)  and 
(43).  In  both  cases,  the  contribution  of  the  cross  terms  is  small. 


29 


5.  CONCLUSION 


The  previous  sections  in  this  report  present  an  analysis  of  the  effects 
of  I  and  Q  Rain  mismatch,  I  and  Q  timlnR  mismatch,  and  timing  jitter  in  the 
sampler  of  pulse-Doppler  radars.  The  simplest  case  is  gain  mismatch  which 
causes  a  false  target  to  appear  at  the  same  range  as  the  true  target,  but 
with  the  opposite  Doppler  shift.  Gain  mismatch  also  causes  an  SNR  loss. 
Sample  timing  mismatch  causes  similar  effects  to  that  of  gain  mismatch,  ex¬ 
cept  that  the  extent  of  the  effects  depends  on  target  range  and  Doppler 
shift.  Timing  jitter  is  somewhat  different  from  the  mismatch  cases;  the  re¬ 
sult  is  a  distortion  in  the  Doppler  processor  output  and  depends  on  target 
range  and  Doppler  shift.  Timing  jitter  does  not  cause  a  false  target  to  ap¬ 
pear  unless  gain  or  timing  mismatch  is  present.  The  distortion  gives  a  sig- 
nal-to-dlstortion  ratio  which,  for  system  clock  frequency  jitter,  is  inver¬ 
sely  proportional  to  the  square  of  the  burst  waveform  length.  In  the  case  of 
independent  random  jitter,  the  distortion  has  a  spectral  peak  at  the  Doppler 
shift  of  the  target. 
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Appendix  A 

In  this  approach,  it  is  shown  that  the  mismatch  term,  6,  can  be 
neglected  in  the  sine  and  cosine  arguments  in  Equation  (16),  and  that 
similarly  the  distortion  terms,  r^S  ,  can  be  neglected  in  the  sine  and 
cosine  arguments  in  Equation  (22).  These  simplifications  make  the  analysis 
in  the  report  considerably  easier. 

Equation  (16)  is  repeated  in  Equation  (A. 1)  for  convenience, 

Xnj  *  cos[2irf(mT+6)+Q)G(t,-f/k+6  ,f  ) 

+  jslnl  2it  f  (mT-6  )+0  ]G ( t '  -f  /k-6  ,f  ) 

+  nj(m)  +  jnQ(m)  .  (A. 1) 

In  this  case,  26  is  a  sampling  time  mismatch  and  it  is  very  reasonable  to 
restrict  the  mismatch  to  be  less  than  one  sampling  period,  26  <  1/kx. 

From  the  definition  of  G(t,f)  in  Equation  (7),  the  real  part  of  the 
signal  portion  of  Equation  (A. 1)  can  be  written  as 

cost  2irf  (mT+6)+Q]  sinlirk(  t'  +6)(t-  1 1'  -f /k+6  I  )  ] 

rrk(t*+6)x  1  1  *  (A.  2) 

In  amining  Equation  (A. 2),  it  should  be  noted  that,  in  pulse-Doppler 

radars,  the  Doppler  shift,  f,  should  be  much  less  than  the  bandwidth,  kx . 

Using  this  relationship  and  the  restrictions  on  t'  and  6,  the  6  component  to 

the  argument  of  the  sine  in  Equation  (A. 2)  is 

tt  k6  (x  —  j  t *  — f  /k-f6  j)  =  itkx6  .  (A.  3) 
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The  6  component  of  the  argument  to  the  cosine  In  Equation  (A. 2),  however,  Is 
2nf  6  .  (A. 4) 

We  can  now  state  the  relationship  between  these  arguments  as 

2n f 6  «  2irkT6  <  u  .  (A. 5) 

Since  the  6  component  of  the  argument  to  the  cosine  Is  much  less  than  that  of 

the  sine  In  Equation  (A. 2)  and  also  much  less  than  it,  It  Is  clear  that 

omitting  2rr  f  6  in  the  cosine  argument  In  Equation  (A. 2)  will  no.  significantly 
affect  the  analysis  of  the  signal.  Similarly,  the  2tt f 6  term  can  be  dropped 
from  the  sine  argument  In  Equation  (A. 1)  giving  the  desired  result. 

The  justification  for  dropping  the  term  itf(rm±6)  from  the  sine  and 
cosine  arguments  in  Equation  (22)  is  essentially  the  same  as  the  discussion 
for  Equation  (16)  given  above.  The  main  difference  Is  that  the  sampling  time 

jitter  plus  mismatch  should  be  less  than  one  sampling  period  so  that 

jr^tfi  j  <  1/k.r  . 
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sampled  inphase  and  quadrature  components  of  the  pulse  compressor  output  must  be  matched  in  order 
to  minimize  distortion  and  noise  at  the  Doppler  processor  output.  In  addition,  the  jitter  in  the  sam¬ 
pling  times  also  causes  distortion  and  noise.  This  report  presents  an  analysis  of  these  effects  and  de¬ 
scribes  the  signal  degradations  due  to  mismatch  and  jitter  in  the  sampler. 
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